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Cotunneling is an important error process in the application of single electron tunneling devices for metro- 
logical and electronic applications. Here we present an experimental investigation of the theory for adiabatic 
enhancement of cotunneling by coherent microwaves. The dependence is investigated as function of temper- 
ature, gate voltage, frequency, and applied microwave power. At low temperatures and applied power levels, 
the results are consistent with theory, using only the unknown damping in the microwave line as a free pa- 
rameter. However, the results indicate that the effects of temperature, frequency and microwave power are not 
independent, contrary to what is suggested by theory. 

PACS numbers: 73.23.Hk, 73.40.Gk, 73.40.Rw 



I. INTRODUTION 

The unique properties of single electron tunneling (SET) 
devices has made them the subject of extensive research. This 
research includes applications such ag current standardsjll ca=. 
pacitance standard^B electrometers,Q quantum computingj^ 
and thermometry .ED The first-order behaviour is adequately 
described by the well-established orthodox theorylM The so- 
called cotunneling, a second-order phenomenon in which one 
electron tunnels through each of the SET junctions at the 
same time, has also received some attention, as this con- 
stitutes an important error process in many SET applica- 
tions. The dependence of the cotunneling current_on volt- 
age and temperat ur a has , b een studied theoreticallyO as well 
as experimentally.ll2llJIHLj Furthermore, Covington et fl/.Lll 
have studied the frequency dependence for 4- and 6-junction 
pumps. This paper investigates experimentally the theoretical 
prediction by Flensberg|i3 that the cotunneling current should 
depend not only on temperature and voltage, but also on the 
amplitude and frequency of an applied oscillating field. 

The paper is organized as follows: The theory for cotunnel- 
ing derived by Flensberg and others is briefly summarized in 
section II. The measurement setup is thoroughly described in 
section |l|. The experimental results on cotunneling are pre- 
sented in sectio n |ry| , and in section ^ the results are discussed. 
Finally, section [Vl[summarizes our conclusions. 



Here, R12 are the tunneling resistances of the left 
and right electrodes, respectively. A* = (e/Ci)(e/2 =F 
CgVg) (mod [e^/Cz]) are the energies to add/remove one 
electron to/from the island, Te is the temperature of the elec- 
tron system, and /?k = h/e^. In the maximum blockade 
state. A* reduce to the charging energy Eq = e^/2Cz, where 
Cz = Ci +C2 + Cg is the total capacitance of the island (the 
ground capacitance is assumed negligible). The validity of 
Eq. |l]has been verified experimentalbdn metallic systems by 
several researchers: Geerlings et ali^ first reported the ob- 
servation of cotunneling current and its scaling with the con- 
ductances and voltage. The scaling of the current, rather than 
the quantitative value, made the results clearly distinguishable 
from thermally enhanced,sequential tunneling. Later, the ex- 
periments by Files et aZ.llil confirmed the temperature and gate, 
bias dependence quantitatively. Furthermore, Pasquier et al\13 
and other groups have investigated cotunneling in 2DEG sys- 
tems. Since it seems well-confirmed by experiments, Eq. |l] 
will be assumed valid here. 

Flensberjli has extended the analysis by Averin & Nazarov 
to the case where a harmonically varying signal Vac cos(27t/f ) 
is applied on top of the DC bias. His result for the cotun- 
neling current, obtained by expansion to the third order in 
the energies, becomes in the adiabatic (low-frequency, low- 
temperature, low-amplitude) limit 



n. THEORY 

Cotunneling (also referred to as macroscopic quantum tun- 
neling of electric charge or q-MQT) is a second-order process 
by which an electron tunnels through the junctions of the SET 
transistor via an intermediate virtual state.LL3 The inelastic 
cotunneling current in a single SET transistor for low temper- 
ature (kBTg <C A*) and bi^ voltage (eVoc ^ A*) was derived 
by Averin and Nazarov as 
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In restating his equations we assume (as is the case for our 
experiments) that the alternating bias is applied to the left lead, 
while the right is kept at a constant potential. 

Eq. ^ is not surprising, as this is what one would ob- 
tain from Eq. |l] by making the substitution Vdc ^ Vdc + 
yAccos(27t/f ) and averaging over time. However, because of 
the power expansion approach, Flensberg was able to derive 
also the next leading order correction to Eq. ^by an expansion 
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FIG. 1 : The orthodox theory (dashed) and cotunnehng theory (soUd) 
predictions of the zero-bias, maximum Coulomb blockade conduc- 
tance for a SET transistor as function of temperature. The simulation 
parameters are the same as those measured for the sample reported 
on. 



to fourth order in the energies: 
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The temperature, ampHtude and frequency are assumed low 
enough that an expansion in powers of these is appropriate. 
Note that in the case of maximum Coulomb blockade, ij:^} 
vanishes, since then A* = Eq and thus the third multiplicative 
term becomes zero. 

The aim of the work on cotunneling presented here is to ex- 
perimentally verify or disprove Eq. We do this by applying 
microwaves to the device and measuring the differential con- 
ductance by a lock-in technique in zero DC bias. This ensures 
that the current predicted by Eq. || does not contribute to the 
result. Measurements to verify Eq. |^ are planned for the near 
future. 

Using the familiar Master Equation approach for the ortho- 
dox theory, we have calculated the curve shown in Fig. |l]for 
a SET transistor with the parameters measured for the sample 
reported on here (JYU NM3). The figure shows the sequen- 
tial tunneling and cotunneling conductances in zero bias and 
maximum Coulomb blockade, as function of temperature. It 
is seen that the conductances have the same order of magni- 
tude, which means that sequential tunneling should be taken 
into account when interpreting the results. 



III. EXPERIMENTAL SETUP 

The samples are fabricated on 500 jum thick oxidized sili- 
con. Gold pads and leads are made by a standard UV lithogra- 
phy technique. Before fabrication, the wafers are cleaned by 



ultrasound. A standard PMMA(-MAA) two-layer resist sys- 
tem is then applied and exposed by a pattern designed for 50 
nm line with and 100 nm overlap. The line width becomes 
a little larger; 70-100 nm due to the exposure equipment and 
development. After exposure, the pattern is developed, and 
cleaned by reactive ion etching. Two layers of aluminium (40 
and 30 nm thick) are evaporated at an angle, with an inter- 
mediate oxidation in 2 mbar O2 for 3 minutes. After the last 
metal evaporation, the sample is oxidized again in 2 mbar O2 
for several minutes. Finally, the result is lifted off in acetone. 

The samples were measured in the KelvinOx™ dilution re- 
frigerator at the Institute of Physics, DTU. The plastic mixing 
chamber has been replaced by a metal one, and the sample 
is mounted on a cold finger extending into a superconduc- 
ting magnet, which, in order to suppress the superconduc- 
tivity of the aluminium, was set to one Tesla in all measure- 
ments presented here. The dilution refrigerator is set up in a 
shielded room, with the mechanical pumps outside the shield- 
ing. The temperature of the mixing chamber could be mea- 
sured with calibrated germanium and ruthenium-oxide ther- 
mometers. The loaded base temperature is 50 mK. 

The samples were biased symmetrically using locally fabri- 
cated low-noise electronics. The electronics is mounted in an 
RF tight, compact metal box for shielding. For bias, an input 
amplifier adds two incoming signals; in this case a DC bias 
and a small low-frequency modulation provided by a lock-in 
amplifier. The added signal is then symmetrized by by an in- 
verting and a non-inverting amplifier. Finally, the signal is 
voltage divided on each side to a suitable bias level Vb, and 
fed to the sample through two large resistors Rq. These were 
set to 10 MQ. in all measurements presented. A FET input, 
low-noise voltage preamplifier with a gain of 1000 was used 
to measure the resulting voltage Vdc across the sample, and 
the voltage signal is further amplified 10 times by an external 
PARI 13 low-noise amplifier. The current through the sam- 
ple is then calculated as (Vb — Vdc)/2/?b- A small-amplitude 
modulation, with frequency 2 Hz, was applied from the out- 
put of a SR850 lock-in amplifier on top of the DC bias to 
measure dynamic resistances. This results in a current mod- 
ulation across the device of about 10 pA. After amplification 
in the FET preamplifier, the voltage response was fed back 
to the lock-in amplifier for measurement. With the measured 
zero-bias conductances, a 10 pA excitation corresponds in the 
worst case (2 /jS differential conductance) to 50 aW heat in- 
put, and the voltage magnitude of the excitation is 5 /jV, cor- 
responding to about 3% of e/Cz- From the measured noise 
level at 2 Hz, combined with the equivalent noise bandwidth 
of the lock-in amplifier, we estimate the RMS noise to be at 
most 2.5% of the signal responses obtained. 

There are 20 DC connections from the room temperature 
electronics to the sample at the base temperature. Starting 
from the 300 K top flange, there are first 20 Thermocoax(S) ca- 
bles (length approx. 1.5 m) connect from room temperature 
to the 1.2 K level, then 20 superconducting wires in a rib- 
bon to the mixing chamber (for thermal isolation), and finally 
20 Thermocoax cables from the mixing chamber to the sam- 
ple (length approx. 25 cm). The 20 Thermocoax cablesand 
the method of wiring minimize cross-talk and external noise 
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input. The cables are thermally anchored at all temperature 
levels. The Thermocoax cables,provide filtering of the room 
temperature thermal radiation;0 for / = 10 GHz (correspond- 
ing to h f /k-Q = 0.5 K) the attenuation of the cable is about 140 
dB/m. The attenuation in dB increases as the square root of 
the frequency. Using a thermally conductive paste, the sample 
is mounted inside a thick-walled copper cavity and contacted 
electrically by "buckling wires" which tread on the sample. 
The 25 cm Thermocoax cables and the sample holder are en- 
closed in a copper shield at the mixing chamber temperature, 
which shields them from the surrounding 4.2 K radiation. 

The individual shielding of the measurement leads also de- 
creased the cross-capacitances considerably, which in turn 
minimizes cross-talk. The total capacitance to cryostat ground 
was measured to be 840 pF per lead atpCjpom temperature, 
which is in agreement with Zorin's result.U It is worth notic- 
ing that the leakage resistance in these Thermocoax cables is 
quite low when the cryostat is at room temperature (down to 
tens of Mi2), but increases to many Gi2 when cooled down. 
We conjecture that this may be due to water vapor being ab- 
sorbed in the MgO insulation powder at room temperature and 
pressure, but being cryopumped and/or frozen out when cool- 
ing the cryostat. This is important because the resistance of 
the SET transistor can be very high in the Coulomb blockade 
state. A matrix connector at room temperature and a "switch- 
board" (20 connectors which can be connected to any of the 
20 slots connecting to the sample) at the mixing chamber fur- 
ther enables us to choose the very best cables for the criti- 
cal connections. Cross-talk can only take place in the room 
temperature electronics, the matrix board, the superconduct- 
ing ribbon, in the "switchboard" at the mixing chamber, and 
in the on-chip wiring. 

To avoid microphonic pick-up from mechanic vibrations in 
the building (pumps etc.), the whole cryostat can be suspended 
using three inflated rubber tubes. Furthermore, lateral vibra- 
tions are damped with a fourth rubber tube concentric with 
the cryostat. This arrangement proved very efficient; compo- 
nents deriving from the asynchronous pumps (around 49 Hz) 
and other components at around 30 Hz vanish upon inflation 
of the rubber tubes. 

The high capacitances in the filters and Thermocoax ca- 
bles force us to use quite low modulation frequencies. For 
our symmetric biasing setup, the cut-off frequency is about 
(Gd + 1/27?b) / (27tQ), where Rb is the value of the the bias 
resistor, Gd is the differential conductance of the device, and 
Q is the total capacitance per lead (the 4-point biasing means 
that there are two leads contributing to the capacitance on each 
side). Because of the high internal cut-off frequency of the 
device, we find it reasonable to assume a simple resistive be- 
haviour for the SET device at these frequencies. In the worst 
case, Gd is smaller than I/IR^, and in this case the cut-off 
will be at l/{4KRBCe). We have used R^ = 10 MQ. in all ex- 
periments reported here, giving a worst-case cut-off at 8 Hz. 
The modulation frequency was chosen to be 2 Hz, from which 
we have calculated the worst-case amplitude error due to the 
filters to be a few percent at the lowest conductances (about 2 

To apply microwaves to the device, there is a separate coax- 



ial connection to the mixing chamber At room tempera- 
ture, a vacuum tight feedtrough connects an SMA connec- 
tor to a standard 0.05" 50 Q. stainless steel coaxial cable in- 
side the cryostat. A DC break is inserted before the vacuum 
feedthrough at room temperature, as well as before the ther- 
mal anchoring at 4.2 K. From this thermal anchor, the signal 
is fed through a short length of Thermocoax cable to another 
thermal anchor at the 1.2 K level. The Thermocoax cable 
serves as a cold attenuator (about 20 dB) for the microwave 
signal as well as unwanted signals such as room temperature 
radiation. The signal is then connected to the mixing cham- 
ber through a 20 m long superconducting 0.05" 50 Q. Nb coax 
cable, which is thermally anchored at the 0.6 K level. The 
superconductivity of the Nb provides a microwave connec- 
tion while retaining thermal isolation of the mixing chamber. 
Before connecting to the mixing chamber, a third DC break 
is inserted in the line. All DC breaks are "inside-outside" 
breaks, meaning that both the inner and outer conductors are 
interrupted. The final microwave connection to the samples is 
provided by one of the 25 cm Thermocoax cables also used 
for the DC connections. This gives further attenuation of the 
signal. 

As a consequence of the distributed attenuation and possi- 
ble resonances in the microwave transmission line and wiring 
to the SET device, the actual power applied to the sample is 
unknown. It would me meaningless to make a throughput 
measurement at room temperature, since the damping char- 
acteristics are different at low temperature, and the coupling 
to the sample is unknown. Therefore, all power values are 
measured at the leveled output of the microwave synthesizer 
The attenuation of the line at each frequency will be a fitting 
parameter in the data. For a given frequency, it is assumed 
that the actual power delivered to the SET device is a fixed 
fraction of the power from the high-frequency source. 

Pressures and temperatures in the circulation system are 
continuously monitored by a computer, which also controls 
the biasing of the sample and reads the results measured by 
the voltage meters etc. through a GPIB bus. The GPIB bus is 
separated electrically from all measurement electronics by an 
optical link, which galvanically separates the computer from 
the measurement, and also enables the computer to be moved 
out of the shielded room if necessary. 



IV. RESULTS 
A. Device Characterization 

For the particular SET transistor reported on here, the 
DC I-Vdc curve (Fig. |]) yield the device parameters 
Rz=Ri+R2 = 95.0kQ. and Ci = Ci + C2 1030 aF. From 
the Vdc-Vg curves (not shown) we get Cq = 0.94 aF. The set 
of Vdc-Vg curves at different current biases also provided ev- 
idence that the device is symmetric, which we will assume in 
the following. The particular sample was chosen from a batch 
for its charging energy, which also means a higher tunneling 
resistance (for fixed RC product). Our simulations show that a 
high charging energy is more important than a low tunneling 
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FIG. 2: Resonances in the high-frequency line. The SET is used 
as a self-detector by biasing in zero bias, maximum blockade and 
applying a fixed —10 dBm signal from the microwave synthesizer. 
The differential conductance is then recorded as a function of the 
frequency. Notice the break in the frequency scale. 
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FIG. 3: Measured I-Vdc curve for sample JYU NM3 SET #2 in 
maximum blockade, and the sequential tunneling prediction. The 
insets show enlargements of the blockade region (upper left) and the 
high- voltage region (lower right). 
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FIG. 4: Temperature estimates from the zero-bias conductance. 
The fit to the sequential tunneling only (squares) provides an upper 
bound, while the fit to sequential tunneling plus cotunneling (points) 
is a lower bound. The electron temperature should approach the mix- 
ing chamber temperature (dashed line) at high temperatures. Also 



shown are the forms 



T(j = 227 mK) and 



Tpjj (dotted lines; = 260 mK for the upper curve, : 
235 mK for the lower curve). 



quencies for the microwave bias, and to determine them we 
used the SET device as a self-detector It was assumed that 
the frequency dependence of the device response was small in 
the relatively narrow frequency range probed. The resonances 
are determined by measuring the dynamic resistance at zero 
DC bias and maximum blockade (using the method described 
above) in small steps. The method is based on the assumption 
of a response both from the sequential tunneling and cotun- 
neling contributions, as neither of these are linear around zero 
bias. The result is shown in Fig. ^ where some clear peaks 
in the dynamic resistance are seen at various frequencies. The 
results shown in this paper are measured at the frequencies 
marked by arrows. 



B. Temperature of the measurements 



resistance in obtaining a high cotunneling/sequential tunnel- 
ing ratio. 

Since we want to first compare with the theory for zero 
current bias and maximum blockade, it is important to make 
sure that this is actually the case in the measurements. To 
do this, we use the following procedure: First, the maximum 
Coulomb blockade is found by sweeping across a peak in the 
Vdc-Vg curve, and fitting to find the exact location of the peak. 
Then, the zero DC bias point was found in a similar manner 
by recording and /-Vdc curve around zero DC bias and fitting 
the dynamic resistance to find the peak. This procedure was 
carried out before each measurement. 

The microwave connection obviously exhibits resonances 
(see section |l|), at which the signal received by the SET is 
higher It is advantageous for us to utilise these resonance fre- 



An important question relating to any measurement at low 
temperature is: What is the correct temperature? First of all, 
the sample may be far away from the cooling source; in our 
case the mixing chamber of the dilution refrigerator This 
means that any heat delivered to the sample by Joule heat- 
ing or radiation may give a higher temperature at the sample 
than at the mixing chamber. Second, there is the question of 
whether the electron system is in thermal equilibrium with the 
phonon systemllj This suggests that it is not correct to rely on 
the temperature measured by the thermometers at the mixing 
chamber, even when these are well-calibrated. 

To get an estimate of the real temperature temperature of 
the electron system in the SET transistor, we use the results 
of the differential conductance measurements at maximum 
blockade and zero bias. Using both the cotunneling theory 
(Eq. ^ and the orthodox theory, the temperature has been 
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fitted to make the differential conductance Gd agree. The 
junction parameters were fixed at the values determined in 
the previous section, leaving as the only free parameter. 
The results of this temperature fitting is shown in Fig. ^ as 
function of the mixing chamber temperature Tpj, (which we, 
following Ref. |l^ assume equal to the phonon temperature, 
hence the symbol). The temperature has been fitted both to the 
sequential tunneling conductance alone (squares), which will 
give an upper bound on the temperature, and to the tunneling 
conductance with cotunneling included (circles), which will 
give a lower bound under the assumption that other effects 
can be neglected. Thus in the last case we have implicitly as- 
sumed that we can use the theory for cotunneling by Averin 
& Nazarov at low temperatures, which has been confirmed by 
several researchers.liJiil Also shown is the "strong coupling" 
curve Te = Tph (dashed line), which can be considered a "hard" 
lower bound. 

It is seen that the fitted electron temperatures indeed sat- 
urate at low mixing chamber|-temperatures, as suggested in 
the work by Wellstood et alx3 However, attempts to make 
the results fit to the 



Tpji form, which applies for 
uniform heating in a thin film, seem to fail for our SET struc- 
ture (To is the temperature at which Tg saturates for Tpj, — > 
0). A much more reasonable fit is provided by the form 
Te = yTo+T^ (solid line). This form respects the lower 
bound as well as the limit at Tph ^ 0. Such a dependence 
is not entirely unreasonable, considering that the heating is 
localized in the SET transistor, and that heat transfer is pro- 
vided only by the leads, which should probably be considered 
between one- and two-dimensional. Indeed the experiments 
of Wellstood et al. on local heating of a SQUID with cooling 
fins gave a T^-"^ form. 
We will adopt the 



- T^^ form as our best esti- 
mate of the temperature in the following. We make no claims 
that this is the accurate temperature, only that it seems to be 
a reasonable approximation which makes a sensible transition 
from the low- temperature to the high-temperature regime. 



C. V^p dependence 

In this section, we present measurements at maximum 
blockade, which means that A* reduce to the charging en- 
ergy Ec- Also, since there is no asymmetry in the device at 
zero bias and maximum blockade (no tunneling direction is 
favored), there should be no frequency dependence in the dif- 
ferential conductance, which is also predicted by theory. The 
absence of frequency dependence enables us to find the rela- 
tive damping of the line at the individual frequencies. Without 
frequency dependence, the zero bias conductance Gd as func- 
tion of the microwave amplitude Vrf should be the same at 
all frequencies. This means that the relative damping can be 
found by minimizing the mean square difference between the 
Gd-Vrf curves at the individual frequencies, using the damp- 
ing as free parameter This procedure provided us with the 
relative damping values 0.0, 3.8, 9.1, 11.5, 14.8 17.3, and 
17.6 dB for the frequencies marked in Fig. ^. As expected 
from Fig. ^ the damping increases with frequency. However, 



CO 



o 



■ 

• » ... 


■ 


■ 


■ 










L(J..^'---.»'^-- 










■\ ■ ■* 

t ^ 


.^^^ ^ 


s..-»-». 




« 

. 






-„-"•' 








♦» ».s — • 
-: 






. ■ 










,1 •' 















20 40 60 80 100 120 140 160 

RF amplitude l/pp [^V] 

FIG. 5: Differential conductance as function of the microwave am- 
plitude at 37.9 dB additional damping. The measured values (points) 
are connected to guide the eye. The lines are the measured zero-bias 
value plus the predicted excess conductance with cotunneling (solid) 
and without cotunneling (dashed). The predictions are at 232, 255, 
275, 301, 336, 379, 428, 478, and 549 mK, respectively (bottom to 
top). 



it should be noted that the optimal values vary slightly with 
temperature, suggesting that there might be a co-dependence 
on temperature and frequency. 

The additional damping (i.e. the damping that should be 
added to the relative damping to get the absolute damping) 
should be the same at all frequencies and temperatures. The 
absolute damping is unfortunately unknown. All we can do is 
to find a damping value that is consistent with theory. In Fig. 
||, the resulting theoretical prediction of the Gd-Vrf curves at 
different temperatures are shown along with the measured val- 
ues, assuming an additional damping of 37.9 dB. It is seen that 
the curves from different frequencies have indeed collapsed 
into one, especially at the lowest temperature. At higher tem- 
peratures, where the change in tunneling is more modest, the 
picture is not as perfect. Also shown is the predicted excess 
conductance with and without cotunneling at the same tem- 
perature added to the zero-power conductance. It is seen that 
the measured conductance fits the predicted one quite well, 
especially at low temperatures. 

The results just presented could lead to the hasty conclu- 
sion that the theory for coherent photon assisted cotunnel- 
ing is coiTect. However, there are several points which cause 
concern. First, the additional damping required to make the 
curves fit changes with temperature, which it should not, un- 
less the change from about 200 mK to about 500 mK really 
increases the damping (which seems very improbable). There 
could be a range of other explanations for this result. For ex- 
ample, the change could be an artefact of the fitting procedure: 
If the excess cotunneling conductance becomes smaller than 
predicted by theory at higher temperatures (which seems rea- 
sonable), this will show up in the fit as a higher damping. A 
related problem is that the quality of the fit should be the same 
at all temperatures if the temperature used for the theoretical 
prediction is correct (i.e. equal to the actual electron tempera- 
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FIG. 6: Differential conductance as function of microwave ampli- 
tude at 36.2 dB additional damping. The data (points) are the same 
as the bottom curves in Fig. |^ except that a lower damping in the mi- 
crowave line has been assumed. The lines are the measured zero-bias 
value plus the predicted excess conductance with cotunneling (solid) 
and without cotunneling (dashed), for a temperature of 257 mK. This 
is the electron temperature fitted from the maximum blockade con- 
ductance (compare the squares in Fig. Note that the horizontal 
scale is different from that in Fig. n. 



6.0 




1.5-1 . 1 . 1 . 1 . 1 . 1 . 1 . 1 . 1 . 

60 80 100 120 140 160 180 200 220 240 

Gate voltage Vq [mV] 

FIG. 7: The A='= dependence (points) and the predicted excess con- 
ductance added to the measured zero-bias conductance, with and 
without cotunneling (solid and dashed lines, respectively). The pre- 
diction with cotunneling is made at 272 mK, which is the assumed 
electron temperature (T^'^ form; line in Fig. m. The sequential tun- 
neling prediction is made at 286 mK; this is the temperature that 
makes the the maximum-blockade point fit (like the squares in Fig. 
^). Also shown is the sequential tunneling part of the prediction with 
cotunneling (dotted line). 



ture). This is because the temperature does not (in the theory) 
contribute to the excess cotunneHng conductance; it only gives 
a constant contribution (compare Eq. However, the fit is 
best at low temperatures. 

One 'sanity check' that should be performed is whether the 
orthodox theory alone can explain the results. Using the value 
36.2 dB for the additional damping gives the result in Fig. ^ 
It is seen that the results cross the prediction for the ortho- 
dox tunneling at this damping value. At higher damping, the 
conductance would be too high at high amplitudes, while a 
lower damping the conductance would be too low at low am- 
plitudes. Thus the results can not be explained by the orthodox 
theory, even if one assumes a higher electron temperature and 
another damping. The electron temperature would also have 
to be substantially higher than at the mixing chamber at all 
temperatures, as seen in Fig. ^ Note that although it would 
seem that we are applying a double standard here, since the 
horizontal scale is different in Fig. ^ and |[ However, the 
orthodox theory should be valid for all voltages and thus am- 
plitudes, while the cotunneling theory only claims validity for 
eyRF<A±. 



D. The A='= dependence 

The dependence of the cotunneling on A* has been inves- 
tigated by measuring the dynamic conductance at zero bias as 
function of Vq. Zero bias was found as described previously, 
and the device was biased within one e-period of the gate bias. 
The result is shown in Fig. ^ along with the predictions of 
the excess conductance with and without cotunneling added 
to the maximum-blockade result. For the prediction including 



cotunneling, the assumed electron temperature is used (solid 
line in Fig. ^; for the sequential tunneling the temperature 
fitting the maximum-blockade point (squares in Fig. As 
expected, the cotunneling prediction fits reasonably well near 
maximum blockade, but diverges quickly when approaching 
minimum blockade, where A* 0. Also shown is the se- 
quential tunneling contribution to the current; it is seen that if 
the temperature used is correct, the sequential tunneling more 
or less 'takes over' as A* approaches zero. The sequential tun- 
neling alone seems to fit surprisingly well. One explanation 
for the good fit is that as_A='= approaches zero, the threshold 
voltage becomes smaller,!!]! and thus the orthodox theory will 
account for most of the current. 



E. Photon assisted cotunneling away from blockade 

The theory derived by FlensbergEl applies for 
eVuc, k^Te, cVrf, hf ^ A*. Thus, when one of A* be- 
comes smaller, the voltage/temperature/amplitude interval 
where this assumption is valid becomes smaller. Measure- 
ments similar to those presented in the previous sections 
have been performed also at CaVc/e = 0.25,0.50,0.75, i.e. 
A^/^c = 0.5,0.0,0.5. Thus the bounds of vaHdity for Vdc, 
and Vrf are diminished at CaVc/e ~ 0.25,0.75 (slopes in the 
Vdc-Vg curve), while the condition is impossible to satisfy at 
CcVa/e ~ 0.50 (minimum blockade), where A* = 0. Another 
way of stating this is that Eq. ^ diverges at A* 0. 

The results are presented in Fig. |[ It is seen that on the 
slopes, the observed conductance falls between the predic- 
tions with and without cotunneling, which is expected since 
the cotunneling should be overestimated at these values of 
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A*. At minimum blockade the conductance approaches the 
prediction from the orthodox theory alone, which is reason- 
able for a SET transistor in the minimum blockade, where the 
sequential tunneling dominates (this is similar to approaching 
the charging voltage at maximum blockade). 



V. DISCUSSION 

The interpretation of the results involves several explicit 
and implicit assumptions. To determine the lowest temper- 
atures, it is assumed that both the orthodox theory for sequen- 
tial tunneliflg and the original cotunneling theory by Averin 
& NazarovO were correct within the given bounds of validity. 
The orthodox theory has been extensively confirmed since the 
first theories were,puhlished, and also the original cotunneling 
has been verified.ll^lU It is further assumed that the coupling 
from the electron system to the mixing chamber through the 
phonon system was strong enough at high temperatures to as- 
sume a convergence towards the mixing chamber temperature 
here. 

The largest problem in the interpretation of the results is 
that both the electron temperature of the device and the damp- 
ing in the microwave line were not known exactly. Since it 
is certainly insufficient to use the mxing chamber tempera- 
ture, which is not in equilibrium with the electron temper- 
ature, it is necessary to rely on existing theories and some 
assumptions to get an estimate of the electron temperature. 
Similarly, it is necessary to rely on the absence of frequency 
dependence in the device to determine the relative damping 
values at the different frequencies. The additional damping 
is then an adjustable parameter Here it is implicitly assumed 
that the damping does not depend on amplitude. Ideally, it 
should also have been checked thoroughly that the differential 
conductance measurement was independent of the modulation 
amplitude. However, this was only done sporadically, i.e. by 
halving and doubling the amplitude and checking that the re- 
sponse followed. 

Another problem is that the temperatures did not range to 
below about 7c/4 (another way of saying this is that the tunnel 
capacitances were somewhat too high). Thus it is difficult to 
say with conviction that the condition feeTe ^ A* is satisfied, 
even at the lowest temperatures. 

Among the other assumptions are that the the SET transistor 
exhibited a resistive behaviour at low frequencies (should be 
reasonable since RzCz corresponds to about 1.7 GHz), that 
the damping in the microwave line is linear, that the leads are 
equal (small correction), and that other effects, such as noise 
and the environment, can be neglected. 

Under the given assumptions, the damping could be chosen 
such that results are consistent with the theory at low tempera- 
tures. Furthermore, it is shown that sequential tunneling alone 
could not explain the results. Thus, the results speak in favor 
of the cotunneling theory. At least we can say that more than 
the orthodox theory is needed to explain the results, and that 
cotunneling is a resonable explanation. However, due to the 
dicrepancies observed, it is clear that a more complete theory 
is needed. 



There is still the unsettling matter of the required additional 
damping given by the fitting procedure is changing with tem- 
perature. A number of explanations can be offered for this. It 
could be the coupling that changes, but it is hard to see how 
such a small change in temperature (from 200 mK to 500 mK) 
would cause a change of 10 dB in the damping. A more rea- 
sonable explanation could be that the enhancement of cotun- 
neling by a coherent source must be temperature dependent 
(i.e. if the real temperature contribution becomes smaller than 
predicted by the theory, this will in the fitting be compensated 
by increasing the damping, making the microwave contribu- 
tion smaller instead), which would require the assumed elec- 
tron temperatures to be approximately correct. For example, 
we may conjecture that at higher temperatures, the limit of 
validity is not far away; i.e. that we should require something 
like (iKhnnf + {eVDcY- + jieVRpf « (A±)2 in Eq. |. An- 
other possibihty is that the temperatures used are too low, with 
the error increasing with temperature. This would make the 
orthodox prediction too high, which would then be compen- 
sated in the fit by a higher damping. However, the orthodox 
theory prediction approaches the conjectured electron temper- 
atures at high temperatures, making the error margin upwards 
smaller. 

Also the relative damping varied with temperature, hint- 
ing that also the frequency might have a different effect on 
the conductance at different temperatures, again contrary to 
the theoretical prediction. Indeed the charging frequency of 
the device is fc = 18.8 GHz, which means that the frequen- 
cies used range from 0.11/c to 0.28/c. Similar arguments as 
above could be applied to this matter. 



VI. CONCLUSION AND OUTLOOK 

The enhancement of cotunneling by a coherent microwave 
source was investigated experimentally and compared with 
theory. The results on cotunneling seem to confirm the theory 
for adiabatic enhancement of cotunneling by a coherent mi- 
crowave source. However, the presence of assumptions and 
adjustable parameters should inspire modesty in the strength 
of the conclusions drawn. Furthermore, the results seem to 
suggest that the effect of microwave amplitude and temper- 
ature are not independent, contrary to what is suggested by 
theory. Thus a more complete theory is needed. Experimen- 
talists would very much appreciate a theory which gives a 
method of actually calculating the cotunneling current at all 
voltages, temperatures etc. (e.g. something like the Master 
Equation approach). In any case, the presence of microwave 
induced cotunneling enhances the importance of taking cotun- 
neling into account in the calculation of SET devices. 

It is vital in future measurements to have samples with 
a high charging energy to make the charging temperature 
higher. Also one could attempt to lower the electron tempera- 
ture, e.g. by using cooling fins on the sample. Cold resistors 
and cross-correlation voltage measurements should be used to 
improve the noise. A new voltage bias electronics, which will 
reduce the noise and eliminate the problem with resistive bi- 
asing of a high-impedance sample, is under construction. The 
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measurements should focus on verifying the photon-assisted 
current in Eq. ^, e.g. by biasing in zero voltage and observing 
the current as function of amplitude and frequency. A care- 
ful determination of the device parameters and damping of 
the microwave line is essential, as there will be adiabatic con- 
tributions from both sequential tunneling and cotunneling. It 
should futher be checked that the differential conductance re- 
sult is idependent of the amplitude of the lock-in modulation. 
Also the measurements should be performed on many more 
samples to rule out sample-specific effects. 
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FIG. 8: Microwave assisted cotunneling away from blockade, using 
the results at the seven frequencies collected by applying their rel- 
ative damping. The curves are for the SET transistor biased at the 
negative slope of the Vbc-^G curve (top), minimum blockade (mid- 
dle) and positive slope (bottom). Note the different vertical scale in 
the middle graph. The temperatures, line types etc. are the same as 
in Fig. H 



